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ABSTRACT 


The study of the dynamics of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) 1 is extended to emphasize the synthesis of control laws for both 
the linearized system as well as the large amplitude slewing maneuvers required 
to rapidly reorient the antenna line of sight. For control of the system 
through small amplitude displacements from the nominal equilibrium position 
LQR tecliniques are used to develop the control laws. Pontryagin's maximum 
principle is applied to minimize the time required for the slewing of a 
general rigid spacecraft system. The minimum slewing time is calculated 
based on a quasi- linearization algorithm for the resulting two point boundary 
value problem. 2 The effect of delay in the control input on the stability of 
a continuously acting controller (designed without considering the delay) 
is studied analytically for a second order plant. System instability can 
result even for delays which are only a small fraction of the natural period 
of motion. 3 
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0 Geoccnter 


FIGURE III-l: 


the 3-D geometry of the scole 
CONF ICURAT ION IN ITS DEFORMED STATE 


. Equations of Motion 
'Derived using a Newton-Euler approach 
. Assumptions 

- Reflector and Shuttle rigid 

- Mast has constant cross-section 

- It is assumed to undergo small elastic deformations only 

- Its modal shapes in orbit are assumed to be the same as 
those of an identical* non-rotating beam. 

Stability Analysis (Rigidized SCOLE) 

A stability analysis of the rigidzed SCOLE was conducted for 
the following configurations: 

Rigid - no offset . Pitch motion decouples from roll 
and yaw. in the linear ranges. System not stable 

b) Rigid - with offset parallel to roll axis. Pitch motion 
still decouples from roll and yaw in the linear range. 
System unstable. 

c) Rigid - With both offsets (parallel to roll and pitch 
axes). The motions in all 3 degrees of freedom are 
coupled. System found to be unstable. 

. Control Laws ■ 

Assumption: All the states of the system are available. 

It was suggested by J.G. Lin that an intuitively appealing 
practical approach to achieve the LCS pointing objective 
is a two-stage procedure, (a) Slew as if riaid then, 

(b) damp-out flexible dynamics. 

.The linear regulator theorv used here to control 
-the linear model of the rigidized SCOLE, 

-The linear model of the actual SCOLE configuration including 
the first fo».r flexible modes of the mast. 

Next 

Preliminary slew maneuvers o f rigidized SCOLE. 
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SCOLE: TRANSIENT RESPONSES 

6.0 *6. INITIAL femnSATIim I» BOLL r AMM EOUCLlBB.' 
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RIGID SCOLE - CONTROL EFFORTS 
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linear model of scole with flexibility 

TlttNS. (ESP. TO AN INITIAL Bfeg. IN ROLL g^Ua. 





LINEAR MODEL OF SCOLE WITH FLEXIBILITY 







SCOLE: TRANSIENT RESPONSES 

6.0 DEB. INI HU. PEBUfiflAIUBi IN PITCH ft 20*1 Biptlu^Aiu ft 









SCOLE: TRANSIENT RESPONSES 

8.0 DEG. INITIAL PERTURBATION IN YAK ftLDM £^aiu^tjK 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 

nut. ibp. n » mm. mnun.n ram ^%>he>)u>Vi'iUIum 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 

raw. op. to a i.oft.oisT.n woe « 





LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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LINEAR MODEL OF 
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LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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MINIMUM TIME ATTITUDE SLEWING 
MANEUVERS OF A RIGID SPACECRAFT 

OBJECTIVE 

DEVELOP COMPUTATIONAL TECHNIQUES TO SLEW 
A GENERAL RIGID SPACECRAFT < INCLUDING 
RIGIDIZED SCOLE ) FROM AN ARBITRARY INITIAL 
ATTITUDE TO A FINAL REQUIRED ATTITUDE 
PRECISELY, AND SATISFYING THE FOLLOWING 
CONDITIONS: 

. IN MINIMUM TIME 

. THE CONTROLS HAVE SATURATION LEVELS 



METHODOLOGY 


THE MAXIMUM PRINCIPLE FROM OPTIMAL CONTROL THEORY IS 
APPLIED TO THE EULER' S DYNAMICAL EQUATIONS AND THE 
QUATERNION KINEMATICAL EQUATIONS OF THE SYSTEM TO 
DERIVE THE NECESSARY CONDITIONS FOR THE CONTROLS. 

THIS LEADS TO THE TWO-POINT BOUNDARY- VALUE PROBLEM. 

AN INTEGRAL OF A QUADRATIC FUNCTION OF THE CONTROLS 
IS USED AS A COST FUNCTION, BUT THE INTEGRATION 
PERIOD OF THIS INTEGRAL, CALLED THE SLEWING TIME, 

IS TO BE CHANGED UNTIL IT REACHES ITS MINIMUM VALUE. 

THE RESULTING TPBVP IS SOLVED BY A QUAS I LI NEARI ZAT I ON 
ALGORITHM < METHOD OF PARTICULAR SOLUTIONS ) . 

EULER' S EIGENAXIS ROTATION THEOREM IS USED TO 
APPROXIMATELY DETERMINE THE INITIAL VALUES OF THE 
COSTATES AND THE SLEWING TIME AS WELL AS THE 
NOMINAL SOLUTIONS WHICH ARE UESD TO START THE 
QUASILINEARIZATION ALGORITHM. 

THE MINIMUM SLEWING TIME IS DETERMINED BY SHORTENING 
THE TOTAL SLEWING TIME UNTIL AT LK**ST ONE OF THE 
CONTROLS BECOMES A BANG-BANG TYPE. 
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CONCLUDING REMARKS 


THE SLEWING MOTION NEED NOT BE RESTRICTED TO A 
SINGLE AXIS MANEUVER. 

. THE GUESSED STARTING VALUE OE THE SLEWING TIME IS 
VERY CLOSE TO THE CONVERGED VALUE FOR THE SCOLE 
EXAMPLES AND SUBROUTINE USED HERE. 

. THE GUESSED INITIAL VALUES OF THE COSTATES ARE 
ADEQUATE FOR THE ALGORITHM TO CONVERGE. 

. THE METHODS USED HERE MAV BE IMPLEMENTED FOR 
PRACTICAL CONTROL SOURCES WHICH MAY HAVE MORE 
CONSTRAINTS . 

AN EXTENSION TO THE MINIMUM TIME SLEWING MOTION 
OF THE SCOLE MODEL CONTAINING BOTH RIGID AND 
FLEXIBLE COMPONENTS IS PLANNED. 
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Appendix - Chapter II 

Stability Analysis of Second Order System 
with Delayed State Feedback 

As a second order differential equation describes the dynamics of 
a single mode of any large space structure, the stability analysis of such 
a system with delayed state feedback is analyzed and the amount of delay 

that can be tolerated by the system without becoming unstable is arrived 
at analytically. 

The differential equation of second order with state feedback can be 
written as : 

*i + 2^-^ * Ui 2 x t - -k^ft-h) - kpijft-h) (1) 

where 

= i^ 1 modal coordinate 
■ i^ 1 natural frequency 
= i^ 1 mode inherent damping ratio 
kj. = rate feedback gain 
kp 3 position feedback gain 
h = time delay 

The feedback gains kp, kp are designed for the required stability and 
transient response specifications without taking the delay into consideration. 

The inherent damping ratio, and the feedback gains, and kp; will 
give rise to five possible combinations as shown in Table 1 and are thus 
analyzed separately for mathematical convenience and easy understanding. 
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Case I: - 0, 1^ - . ... 


- 0 and k > 0 
r 

The differential equation of the system can be written as: 
*i + “i x i ■ - k A(t-h) 


C3) 


Case 


I 

- 0 

> 0 

=* 0 

II 

> 0 

> 0 

■ 0 

III 

- 0 

> 0 

> 0 

IV 

> 0 

- 0 

1 0 

V 

> 0 

> 0 

f 0 


Note: The remaining three combinations are 
neither feasible nor of interest. 

Table 1 : Feasible Combinations of?!, v k 
for Stability Malysis 1 * * P 

and the corresponding characterstic equation is given by: 

s 2 + u? + Zc^se’ 8 * 1 = o 


k r = 2? i"i* 


C4a) 


where 

The value of h for which the roots of equation (3) cross the imaginary 
axis can be evaluated by substituting s = j M . 

Thus ^ 2 

“i + ^ 25 i w i“Sinwh + 2c.u.ucosuh = 0 .... 

i (4bJ 

For equation (4b) to be satisfied 

sintuh = 0 

and 


2 2 _ 

“i ' u + cosdih = 0 


( 5 ) 
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Thus 


uh = V 2 


and 


h = 


v 2 


«i + • / ”l+; i 2 ] 

Case II: 5! >0, k r = 2?^ and kp = 0 

The characteristic equation of the system described by equation (1) 
is given by 

Ca.?- U 2 +2c i< u i a)Sinmh) + j ^.u^C^coScoh) = 0 

Thus cosuh * -c !/c . 

r i 

and h » 


(6) 


(7) 


cos 


" 1 C i: i/ g i 


“i [ / 'i-'i 2 + ' 3 

For the case where qsq thesystem will always be stable since no value 
of h exists for which the roots of (7) cross the imaginary axis. A plot 
of Wjh versus q for various values of q is shown in Figure 2.1. 


( 8 ) 


Case III: 


C i 


°> kp 


k-r > 0 


The characteristic equation is given by 
.-sh.,_ ,-sh 


or 

Thus 


sZ+ToZ+kySe'^+kpe" 514 . o 

(wl-w^k^inuh+kpcosuh) + j ( u)k r cos<uh- kp S inuh) = 0 
tan tuh * M kr 

^ ^ = 7 [( 2 u?+k 2 ) + 7[j^ + 4<u2k 2 +4k^]] 

Plots of ha>i versus k^ for various values of kp/ w 2 are shown in 
Figure 2.2. It can be seen here that these are many combinations of kp and 
k r for which the roots of Eq. (10) can cross the imaginary axis - i.e. 
value of ha^ which leads to instability. 


(9) 

( 10 ) 


( 11 ) 
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Case IV: 


Thus 


z i >0 > k r 3 0 , kp f o 
We characteristic equation is given by 
(“i 2 '“S c °S“h) * j(2 Ci ' v -k o 

?r ... V 


sinwh = 




“ 2 ' “i^ 1 - 2 ^ * “i 2/ ta-25i 2 ) 2 ♦ V“i 2 ) 2 ] (14) 

Plots of K versus y „2 for various values rf 
Figure 2.3 

Case V: C-* > 0 k >. n u j 

l u » K r > 0, kp f o 

The characteristic equation is given by 
(to i 2 -w 2 +u,k r sinuh+kpcosuh) 

+ j(2s! w i<*i + wk r cosuh- k p s inwh) = o 

By equating the imaginary Dart tn 7*™ , 

senary part to zero, «h can be evaluated as 


' " h * sin" 1 ( ^ 2g i V . , cok- 

^ kr> 2 +, Z 2 " - tan ~ 1 ( — _L ) 

K P -ho k r . kp ' 


after substituting wh in the real part of equation (Is, 

the following equation in the single unknown variable e’can be obtained 

2 2 i 

“i ^Vin (sin'V-tan^C Si 33 

+ k p cos (sin’ 1 y-tan" 1 ( Si 33 = 0 

kp C17' 
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Using equations (17) and (16), the limiting value for given values of C .' 
k r* ^ 311(1 “i can ^ determined. As the equation (17) is nonlinear, 
numerical procedures may have to be used and thus the generalized plots 
similar to the other cases may be obtained; 


\7G 


2-14 



ho) 


i 


Z f 3 inherent damping 
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